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We show that oblique propagation of electrons in crystals of Ge and Si, where the electron velocity does 
not follow the electric field even on average, can be explained using standard anisotropic theory for indirect 
gap semiconductors. These effects are pronounced at temperatures below ~1 K and for electric fields below 
~5 V/cm because inter- valley transitions are energetically suppressed forcing electrons to remain in the same 
band valley throughout their motion and the valleys to separate in position space. To model, we start with an 
isotropic approximation which incorporates the average properties of the crystals with one phonon mode, and 
include the ellipsoidal electron valleys by transforming into a momentum space where constant energy surfaces 
are spheres. We include comparisons of simulated versus measured drift velocities for holes and electrons, 
and explain the large discrepancy between electrons and holes for shared events in adjacent electrodes. 



When low energy recoils occur in cryogenic Ge or Si 
crystal radiation detectors, such as those used for the 
CDMS (Cryogenic Dark Matter Search) experiments,^ 
energy is contained within a few hundred /im as a cloud 
of energetic electrons rapidly relax to very low energy 
initially by exciting more electrons and then by emit- 
ting optical phonons. This cloud has an equal number 
of electrons and holes which begin to separate and move 
in opposite directions in the applied electric field while 
emitting low energy acoustic phononsP In the CDMS 
Ge detectors electrons behave very differently from holes. 
For example, with inner and guard ring charge electrodes 
we find that when electrons are pulled towards the elec- 
trodes the fraction of shared events is ten times greater 
than when holes are pulled towards the same electrodes, 
suggesting ten times larger lateral straggle for electrons P 

The existence of these anisotropic effects have been 
known since the 1950s, but in classical calculations^^ 
only approximations of this physics could be incorpo- 
rated. The computing power available today allows more 
rigorous computation taking into account energy momen- 
tum conservation, the full shape and anisotropy of the 
carrier and phonon dispersion relations, and the vari ous 
scattering processes (e.g. intra- valley and inter- valley) 

However, it is important to develop simpler models 
to gain physical insight of the important processes and 
to validate the more ambitious transport calculations. 
Such simplifications may also speed up Monte Carlo de- 
tector simulations. In this spirit, we present a sim- 
ple model where for holes the anisotropic nature of the 
crystal is taken into account by considering an isotropic 
medium with the average properties of the anisotropic 
crystalsPand additionally for electrons we transform the 
anisotropic momentum space to an isotropic spaceP We 
show that this model successfully reproduces the drift ve- 
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locities of the holes and the electrons using one adjustable 
parameter (deformation potential S) for each. Both Ge 
and Si are indirect gap semiconductors with electron val- 
leys away from the momentum space origin, and the hole 
valence band with its minimum at the origin of momen- 
tum space. For holes, the model assumes a T spherical 
valley, and for electrons the model includes the four L 
elliptical valleys along [111] axes for Ge or the six A val- 
leys along [100] axes for Si by transforming each into 
a momentum space whe re t he constant energy elliptical 
surfaces become spheres Then transforming back to 
real space we obtain the lowest order anisotropic behav- 
ior for the electrons. These anisotropic effects^ dominate 
at temperatures below ~1 K and for electric fields below 
^5 V/cm because inter-valley transitions are energeti- 
cally suppressed forcing the electrons to remain in the 
same valley throughout their motion. 

Starting with an isotropic model for holes, we consider 
the case of an incident hole with energy eg and wave vec- 
tor k scattering off of the lattice and emitting a phonon 
of energy huj and momentum q, and the hole has final 
state energy eg, and wave vector kl . Below 1 K thermal 
phonons are suppressed so that only spontaneous phonon 
emission occurs. Energy and momentum conservation re- 
quires eg — eg, = hw and k — k' — q + K where the recip- 
rocal lattice vector K = at our low temperature and 
low electric field. Since only longitudinal phonon modes 
couple in an isotropic medium, we assume one mode with 
ui/q = sl and isotropic speed of sound sl- We assume 
a parabolic band minimum, so that eg = h 2 k 2 / '(2m^) 
where mj, is the hole effective mass. In Table [TJ we col- 
lect the constants needed for these calculations. 

From momentum conservation and the law of cosines 
k' 2 = k 2 + q 2 — 2kq cos 8 with the angle between k and q, 
and from energy conservation k 2 — k' 2 = (2mi l /h 2 )?no — 
(2rrihs/h)q = 2k s q where k s — m^s/h. Combining we 
obtain q = 2(fccos6' — k s ), which has solutions only if 
k > k s or the hole velocity v is greater than the speed of 
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TABLE I. Physical constants for Si and Ge crystals. The 
isotropic hole effective mass rrih, and the anisotropic elec- 
tron effective masses my and m± are || and _L, respectively, 
to the conduction valley axes, and conductivity effective mass 
3/m c = 1/mii +2/m±. The incident energy per final electron- 
hole pair is e e h, sl the speed of sound, and lo = nh 4 p/ (2m 3 H 2 ) 
is the characteristic range for carrier scattering where Si 
(frorrP') or Sut (fit to dateP) is the deformation potential. 
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sound s, analogous to Cherenkov radiation of photons. 

For simulations,^ we use Fermi's golden rule to com- 
pute the azimuthally symmetric differential phonon emis- 
sion rate P(k, 9) sm9d8 for angles between 9 and 9 + d9 
relative to the incident hole direction 



P(k,9) sin 9d9 = 
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where < 9 < cos" 1 (k s /k) < 7r/2, and = 
irh 4 p I '(2m|H 2 ) is the characteristic scattering range 5 
with deformation potential S. Integrating over 9 we ob- 
tain the total isotropic hole-phonon scattering rate 
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for k > k* 
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Finally, we compute the hole scattering angle and obtain 

k 2 -2kJkcos9 - k s ) -2(fccos6> - k s ) 2 
cos cj) = V s > y tL_ ( 3 ) 

k^/k? - 4k s (k cos 9 - k s ) 

where the angle ranges from < 4> < it for k > 2k s and 
from < <f) < tt/2 for k s < k < 2k s . 

Using these equations for hole propagation in the 
isotropic medium, we have performed simulations at var- 
ious applied voltages. In Fig. [T] we show 20 hole trajecto- 
ries each with 1000 scattering events. The spatial disper- 
sion about the mean displacement is a diffusive process, 
and increases with the square root of distance or time. 
In Table [TT] we show analytic versions of many of the pa- 
rameters of interest in the limit k » k s , including the 
lateral dispersion, average phonon energy, average hole 
energy, and hole drift velocity. Finally, we compare with 
the Ge drift velocity measurements made recently with 
CDMS detectors^ and we fit the simulation to the data 
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FIG. 1. Simulation of hole propagation through Ge cyrstal at 
zero temperature with a 1 V/cm applied electric field. Shown 
are the hole trajectories (green), the emitted phonon posi- 
tions (blue), and the final hole positions (red). On average 
1000 scatters took 0.59 fis, traversed 10.8 mm and transverse 
straggle standard deviation of 0.71 mm. 
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FIG. 2. Comparison of the measured and simulated drift 
velocities as a function of electric field along [100] axis for 
ultrapure (N a = 1.83 x 10 11 cm" 3 ) Ge crystals at 31 mK. 8 - 
One parameter fits to Ref 8 data were obtained by setting 
the deformation potential for holes to H/ it = 3.4 eV (green) 
and for electrons to H/it = 11.0 e V (r ed). Also shown are 
published theory curves for electrons^' and for holesP 



with = 3.4 eV. Note from Fig. [2] that the shape of 
the hole drift velocity versus electric field is reproduced. 

Next, we model the electrons near the minimum of 
their anisotropic valleys. When an electric field is applied 
the rate of change in momentum, eE — dp/dt = hdk/dt, 
is along the direction of the electric field, but because the 
effective mass is a highly anisotropic tensor, p = rhv and 
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TABLE II. Drift velocity parameters for Si and Ge, in terms of 
a unit-less a(E) = (eElo) / (2e 3 ) which is proportional to the 
applied electric field E and e s = |m c s 2 . For isotropic models, 
in the limit of k >> k s , we compute the drift velocity Vd, 
average carrier velocity v, average phonon energy hui, average 
carrier energy e, mean free path A and component parallel to 
field, and standard deviation of straggle a±. The values in 
last four columns are for E = 1 V/cm and x = 1 cm. 
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FIG. 3. Simulation of electron propagation through Ge 
cyrstal at zero temperature with a 1 V/cm applied electric 
field along the -z axis (x, y, & z aligned with [100] crystal 
axes). Shown are the electron trajectories (green) and the 
final positions (red). On average 1000 scatters took 0.28 fis, 
traversed 7.4 mm along the z-axis. There is large anisotropy 
with four elliptical cross section beams with the plane of each 
near perpendicular to one [111] axis (each beam is 33.2 deg 
from z axis). The ratio of elliptical spreading is ~ \ fmvJrnZ. 

(m) _1 ei? = dv/dt, the rate of change in group velocity 
is in general far from the electric field direction. For a 
coordinate system with the z-axis aligned with the prin- 
ciple axis of the valley and origin at the valley minimum, 
we have eEi/rrii — dvi/dt and the energy is given by 
= h 2 kl/(2m 1 _)+h 2 kl/{2m 1 _) + h 2 kl/{2m\\) which are 
azimuthally symmetric ellipses with principal axes ratio 
of y/m\\/m± = 4.41 for Ge. 



For each valley, we transform into a momentum space 
where the equal energy surfaces are spheres with the 
same energyP Thus k\ = Wmi/m c k* and we derive 
Ei = ^/m l /m c E*, v t = v*/y/m.i/m c , Xi = x*/y/mi/m c , 
where 3/m c = 1/tom +2/m±, and u> = lo* (since eg— eg/ — 
e *k ~ € *k'^' ^ote that a uniform field in real space trans- 
forms to a uniform field in the starred space but mag- 
nitude and direction change. The only remaining issue 
is the speed of sound, which if isotropic in real space is 
anisotropic in the starred space. We assume the speed of 
sound is isotropic in starred space, using the value from 
real space. Within this framework, we perform the entire 
simulation for electron propagation in starred space with 
starred electric fields, then we transform back into real 
space rotating the z-axis back to the four [111] valleys. 
The simulation equations are identical to Eqs [T]|3] where 
all parameters are replaced with their starred partners. 

A typical result is shown in Fig. [3] where we have pop- 
ulated the four Ge valleys with 80 electrons total and 
follow their trajectories in real space with the electric 
field along the [100] crystal axis (-z axis). After the first 
few mm's, the electrons in the four momentum space val- 
leys become separated in real space. The spread for each 
valley around each of the four drift velocity directions 
has a similar average value to the isotropic hole spread, 
but the shape is elliptical with the ratio ~ \/m\\/m±. 
As shown in Fig. [2j the simulations are compared with 
the drift velocity measurements for electrons^ and fit the 
data with Hgt = 11.0 eV. Again the shape of the electron 
drift velocity versus electric field is reproduced. 

Oblique propagation of electrons is absent in [100] Si 
crystals because the six Si valleys all have one principle 
axes aligned with the electric field, but would appear for 
[111] Si crystals. For Ge these effects are reduced in [111] 
crystals, but remain because only one of the valleys has 
principle axes aligned with the electric field. 

In conclusion, we now understand the large spread- 
ing in the propagation of electrons versus holes in Ge 
crystals at temperatures below ~1 K and electric fields 
below ^5 V/cm as due to the freeze out of inter- valley 
scattering in the highly anisotropic electron bandsPNote 
that when transitions between valleys turn on, the sum 
of all valleys becomes isotropic. We thank A. Bronia- 
towski for valuable discussions on oblique propagation. 
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